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Evidence of massive neutrinos in the n + — >• e + v decay spectrum was sought with the background 
H + — > e + decay chain highly suppressed. Upper limits (90 % C.L.) on the neutrino mixing 



ty — >• \i — > e 
matrix element It/, 



in the neutrino mass region 60-129 MeV/c were set at the level of 10 



PACS numbers: 13.20.Cz, 14.60.St, 12.15.Ff 



INTRODUCTION 

A natural extension of the Standard Model (SM) in- 
corporating neutrino mass and possibly explaining the 
origin of dark matter involves the inclusion of sterile neu- 
trinos mixing with the ordinary neutrinos Ijj . The weak 
eigenstates v Xk of such neutrinos are related to the mass 
eigenstates Vi by a unitary matrix, vg = Y^^UuVi, where 
I = e, jH, r, xii X2---Xk- An example of a sterile neutrino 
model is the Neutrino Minimal Standard Model that adds 
to the SM three massive gauge-singlet fermions (sterile 
neutrinos) Q. In the context of this model, a search for 
extra peaks in the tt + — > e + v decay spectrum is sensi- 
tive to sterile neutrinos depending on the mass hierarchy 
structure and choice of parameters Q. 

The decay tt+ -> e+v (E e + = 69.8 MeV) with a 
branching ratio of R = (1.230 ± 0.004) x 10" 4 0-Q is 
helicity-suppressed by (m e /m M ) 2 in the SM. The relax- 
ation of this condition for massive neutrinos facilitates 
the search for extra peaks in the lower positron energy 
region. Previous results Q at the level of \U e i\ 2 < 10~ 7 
in the neutrino mass region of 70-130 MeV/c 2 were lim- 
ited by the presence of unsuppressed [i + — > e + vv decay 
background (E e + — 0.5 — 52.8 MeV) originating from 
decay-in-flight of pions 0]. The TRIUMF PIENU exper- 
iment [8] aiming at a more precise measurement of the 
branching ratio R was designed to further reduce this 
background. 

In this paper, we present results of the search for low- 



energy peaks in the background-suppressed spectrum of 
7T+ — > c+u decays at rest. 



EXPERIMENT 

The extension of the TRIUMF M13 beam line H, 
which suppressed positrons in the beam to < 2 % of 
pions, delivered a 75±1 MeV/c pion beam to the exper- 
iment. Figure Q] shows a schematic view of the detector. 
The 7r + beam was degraded by two thin plastic scintilla- 
tors Bl and B2 (6 mm and 3 mm thick, respectively) and 
stopped in an 8-mm thick active target (B3) at a rate of 

5 X 10 4 7T+/S. 

In order to obtain a large solid angle for e + detection 
from the decay n + — > e + v while reducing variations in 
material along the positron path, the primary calorime- 
ter was placed on the beam axis. It consisted of two 
plastic scintillators (3-mm thick Tl and 6-mm thick T2) 
for positron identification and a 48-cm (dia.) x 48-cm 
(length) single-crystal NaI(T£) detector [10( viewed by 
19 photo-tubes. Two rings of 97 pure Csl crystals [lH (9 
radiation length radially) surrounded the Nal(T^) crys- 
tal to reduce leakage of the electromagnetic shower pro- 
duced by a positron. The observed energy response of 
the Nal(T^) detector to positrons from tt + — > 0+^(7) de- 
cays included a peak at 70 MeV and a low-energy tail 
distribution. 
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FIG. 1: PIENU Detector (see text). 



Pion tracking was provided by wire chambers (WC1 
and WC2) each with three-planes (oriented at 0°, ±60°) 
at the exit of the beam line, and four planes of 0.3-mm 
thick X and Y Si-strip counters S1X/Y and S2X/Y lo- 
cated immediately upstream of the B3 counter. Positron 
tracking came from one set of X and Y Si-strip coun- 
ters (S3X and S3Y) immediately downstream of the B3 
counter, and three layers of wire chambers (WC3) in front 
of the Nal(T^) crystal. 

A coincidence of beam counters, Bl, B2 and B3, with 
a high Bl threshold defined the pion signal, and a coinci- 
dence of positron counters, Tl and T2, defined the decay- 
positron signal. A coincidence of pion and positron sig- 
nals with a time window of -300 ns to 500 ns was the basis 
of the main trigger logic. This was prescaled by a fac- 
tor of 16 to form an unbiased trigger (Prescaled trigger), 
and the early time window (2-40 ns) provided another 
trigger (Early trigger) that contained most of ir + — > e + v 
decays. Beam positrons for detector calibration were ac- 
cumulated by a separate trigger that required a low pulse 
height in the Bl counter. A typical trigger rate was 600 
Hz, of which 240 Hz was for Prescaled, 160 Hz for Early, 
and 5 Hz for Beam positron. 

The pulse shapes from the plastic counters, the Csl 
crystals and Si-strip detectors, and the NaI(T€) crystal 
were digitized at 500 MHz, 60 MHz, and 30 MHz, 
respectively. 



ANALYSIS 



Event selection 



The charge and time of each pulse were first extracted 
from the waveform. Energy calibration of each detector 
was based on the energy loss of minimum ionizing par- 
ticles (75 MeV/c positrons in the beam) for plastic and 
Si-strip counters, and 70 MeV positrons from 7r + — > e + v 
decays for the Nal(T^) signals. Gain instability was cor- 
rected on a run-by-run basis (typically, every ten min- 
utes). 

Events originating from stopped pions were selected 
based on their energy losses in Bl, B2, S1X/Y and 
S2X/Y, and the time-of-flight with respect to the 
primary-proton beam burst (or the phase of the cyclotron 
Radio Frequency) . Events caused by beam positrons and 
muons were negligible after the cuts. Any events with 
extra activity in the beam counters including WC1 and 
WC2 were rejected. About 40 % of events survived the 
cuts. 

An energy-loss cut based on the minimum energy- 
loss in Tl, T2, S3X and S3Y was used to select decay- 
positrons. The cut reduced beam muons directly hitting 
those counters and protons coming from pion-nucleus re- 
actions to a negligible level. A fiducial cut, requiring the 
hit position at WC3 to be within 8 cm from the beam 
axis, was imposed to reduce the low energy tail of the 
7r + — > e + v peak due to shower leakage from the NaI(Ti?) 
crystal. The low energy tail was further reduced by re- 
jecting events with energy above 6 MeV in the Csl crys- 
tals. 
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FIG. 2: Total energy in the beam counters. The peaks at 
15.6 and 18.5 MeV are from 7r + — > e + u and n + — > n + — > 
e + decays, respectively. The selected region is between the 



vertical dashed lines. 



In order to enhance the small 



nal, the 7r + 



c + v decay sig- 
c+ (ir + — > [i + v decay followed 
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by ^ 



vv decay) background was suppressed using 



timing cuts to take advantage of the lifetime difference 
between pions and muons (tv = 26ns and r M = 2197ns). 
The time window selected was 2-33 ns following the pion 
stop time. Since the decay 7r + — > c + v involves only two 
charged particles while the decay tt + — > fi + — > e + has 
three charged particles with an extra kinetic energy of 
4.1 MeV deposit from the ir + —> [i + v decay in B3, pulse 
shape discrimination based on the likelihood for two and 
three pulses, and the total energy in the beam counters 
were also very effective in 7r + — > /i + — > c + background 
suppression. In Figure [5J a spectrum of total energies in 
Bl, B2, B3, S1X/Y and S2X/Y integrated over a time 
window of 100 ns is shown. The peaks at 15.6 and 18.5 



MeV are from 



c+ decays, 



respectively. (The energy separation between the two 
peaks is smaller than 4.1 MeV because of saturation ef- 
fects in the plastic scintillator.) The vertical dashed lines 
indicate the cut positions. The ratio of the low energy 
events (E e + < 54 MeV, including the ir + — > e + v tail) 
and the tt + — > e + v peak (E e + > 54 MeV) was 0.2, consis- 
tent with that obtained in the previous TRIUMF exper- 
iment [7] (before the optimization process described be- 
low). At this stage, the major low-energy background in 
the background-suppressed spectrum came from decay- 
in-flight (DIF) of pions near the B3 counter, in which 
the muon from the n + —> fi + is decay stopped in B3 and 
deposited the same kinetic energy as the initial pion jij. 

The tracking detectors, S1X/Y and S2X/Y, allowed 
detection of a kink in the pion track when DIF happened 
upstream of the B3 counter. For the remaining events, 
the "pion" direction near B3 with respect to the beam 
direction is plotted in FigJ3]for the regions E e + > 54 MeV 
(mostly 7T+ ~> e + v) and E e + < 32 MeV (mostly pion DIF 
events with a negligible tt + — > c + v tail contribution) by 
thin and thick lines, respectively. The background was 
suppressed by another factor of two by requiring the kink 
angle to be < 14°. 
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FIG. 3: Kink angles for events E > 54 MeV (thin) and E < 32 
MeV (thick) . The vertical dotted line at 14° indicates the cut 
position. 



Also, consistency tests of the "pion" and positron 
tracks based on the closest approach of the two tracks 
in B3 provided an extra handle for suppressing the pion 
DIF events. An additional suppression factor of three 
was obtained with a loss of statistics of 30 %. 

The cuts were optimized by minimizing the value 

S = V V <5 4Mcv/-^'>54MeV, where jV <5 4McV and V >5 4M e V 

are the numbers of events below and above 54 MeV in the 
positron energy spectrum, respectively. The final back- 
ground (including the low-energy tail of the tt + — > e + v 
peak) to peak ratio was iV<54Mev/A r >54MeV = 0.068 with 
V>54McV = 4.8 x 10 5 . The positron energy spectrum is 
shown as "No cut" in Fig. [4j 
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FIG. 4: Positron spectra with (35°) and without (No cut) 
angle cuts together with the positron beam spectrum shifted 
by 2.5 MeV (shaded). 



Spectra with angle cut 

There was a strong correlation between the emission 
angles of the decay positrons and the amount of the low 
energy tail due to shower leakage. As shown in FigJH 
the peak resolution was also improved by 10 % with the 
positron angle cut at 35° with respect to the beam axis 
(a = 0.8 % (RMS) at 70 MeV with the 35° cut). This 
cut also contributed to a better peak-to-background ra- 
tio. However, tighter cuts on the positron emission angle 
worsened the S value. Using spectra with the angle cut 
was effective only above 47 MeV where the impact of the 
resolution was higher. 

There was a shoulder in the positron spectrum 
approximately at 60 MeV that was enhanced with the 
angle cut. Similar structures observed in the spectrum 
of the 75-MeV/c beam positron shown by the shaded 
histogram in Fig|4] were discussed in Ref . fl2| . (To 
compensate for the difference in the initial energies and 
the additional energy losses of the beam positrons in 
Bl, B2 and B3, the beam-positron spectrum was shifted 
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by 2.5 MeV to line up with the 7r + — > c + v peak.) The 
bumps (or shoulders) at 60 MeV and 53 MeV in the 
beam positron spectrum correspond to the primary peak 
energy minus one and two neutron separation energies in 
127 I, respectively; the loss of the energy observed in the 
NaI(T£) crystal is due to low-energy neutrons produced 
in photo-nuclear reactions escaping from the crystal. 
Since this 60 MeV shoulder in the suppressed spectrum 
is consistent with the response function of the Nal(TY) 
crystal 
study. 



12j . it was treated as background in the present 
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FIG. 5: Normalized amplitudes (%) of potential peaks for 
the no-cut spectrum (open circles) and for the 35° spectrum 
(closed circles). 



from the late time window did not include this compo- 
nent. The muon-DIF spectrum was obtained by applying 
a Lorentz transformation for the muon kinetic energies 
3.3-4.1 MeV to the n + —> fi + —> c + spectrum with ac- 
ceptance corrections (described in the next section) and 
muon polarization effects. The resulting muon-DIF spec- 
trum had a broad bump around 30-40 MeV, extending 
near the 7r + — > e + i> peak. The amplitude of this compo- 
nent was a free parameter of the fit. 

In order to accommodate a slowly changing spectrum 
mismatch and unspecified background, the amplitude 
and the decay constant of an exponential function and 
an additional constant term were free parameters of the 
fit. 

The template peak spectrum for each peak energy was 
obtained by applying the same cuts used in the data anal- 
ysis to the positrons which were generated by a Monte 
Carlo (MC) simulation in the B3 counter with the ob- 
served pion stopping distribution. Measurements using 
positron beams at various entrance angles and energies 
into the Nal(T^) crystal confirmed the validity of the MC 
line shapes including the effects of Csl vetoing. Agree- 
ments between the data and MC in the peak shape were 
within 10 %. 

Figure [5] shows normalized amplitudes of fitted peaks 
for the spectra without the angle cut (open circles) and 
with the 35° angle cut (closed circles). The x 2 /DOF's 
were 0.97 without the angle cut and 1.00 with the angle 
cut for the fits without the extra peak. 



Fitting 

The search for extra peaks was conducted with 0.5 
MeV steps in the positron energy regions, E e + — 10 — 47 
MeV for the spectrum without the angle cut and E e + = 
47 — 60 MeV with the 35° cut to take advantage of im- 
proved energy resolution obtained using the cut. The 
entire spectrum (9-50 MeV for the no-cut data and 9-62 
MeV for the data with the 35° cut) was fitted to a back- 
ground function (described below) plus a possible peak. 

In order to minimize the effects of the bumps in the 
response function, the beam positron spectrum was sub- 
tracted from the spectra before the fitting search. The 
normalization was done in such a way that the 60 MeV 
peak amplitude in the fit of the spectrum with the 35° 
cut was zero. 

The amplitude of the tt + — >• /i + —> e + background, 
for which the spectrum was obtained with a late time 
window (150-500 ns), was a free parameter in the fit. 

The positrons from DIF muons emitted promptly fol- 
lowing tt + — > [i + v decays (about one third of the total 
7r + — s- [i + — > e + background) decreased with the pion 
decay time, and the ir + — » fi + — > c + spectrum obtained 



Acceptance 

Since the low-energy peak amplitudes obtained were 
normalized to that of the 70 MeV peak, most acceptance 
effects canceled to first order, especially those related 
to the pion definition cuts. There were, however, some 
energy-dependent effects in the cuts to be corrected for. 
The acceptances were estimated based on MC calcula- 
tions. The consistency was tested to be within 3 % by 
comparing the MC and experimental tt + — > fi + — > e + 
spectra in the 10-50 MeV region with and without the 
background suppression cuts. 

The fiducial cut increased the relative acceptance for 
10 MeV positrons by 5 % (for no angle cut) with respect 
to 70 MeV positrons due to multiple scattering effects. 
Energy leakage into the Csl crystals for higher positron 
energy resulting in rejection of events also increased the 
relative acceptance of 10 MeV positrons by 15 %. 

Because of larger scattering cross sections at lower 
energy, even within a small path length in B3 (3 mm in 
depth), lower energy positrons tended to have a larger 
total energy deposit in the target, thus lower energy 
events looked more like tt + — > /i + — > c + decays, causing 
a 10 % loss (no angle-cut data) in efficiency due to the 
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total energy cut. The largest energy-dependent effect 
was in the vertex consistency requirement for pion and 
positron tracks, which reduced the acceptance of low 
energy positrons by 60 %. The combined acceptances 
for 10 MeV positrons with respect to 70 MeV positrons 
were 45 % (35° data) and 42 % (no cut). 



CONCLUSIONS 



The present experiment improved the upper limits on 
the neutrino mixing matrix element | U e i | 2 by a factor of 
up to four in the mass region 68-129 McV/c 2 . 



RESULTS 



No significant peaks above statistical fluctuations were 
observed. After correcting for the acceptance and the 
helicity-suppression and phase-space terms, the ampli- 
tudes and associated errors were converted to 90 % C.L. 
upper limits on |J7 e i| 2 , assuming a Gaussian probability 
distribution with a constraint that the physical region 
of a peak area be positive. Figure [6] shows the combined 
results for the fits with the 35° angle cut (below 80 
MeV/c 2 in neutrino mass), and without the angle cut 
(above 80 MeV/c 2 ). The region below 60 MeV/c 2 
(E e + > 57 MeV) was excluded in the plot because of 
the strong bias caused by the background subtraction 
procedure. For comparison, the 90 % C.L. upper lim- 
its obtained in Ref. [7J are also plotted by a dashed curve. 
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FIG. 6: Combined 90 % C.L. upper limits obtained from the 
35° spectrum (circles) and no-cut spectrum (triangles) to- 
gether with the previous limits (dashed line) [3]. 
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